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Abstract 



. Recent cosmological data for very large distances challenge the validity of the stan- 

\^ I dard cosmological model. Motivated by the observed spatial flatness the accelerating 

■ expansion and the various anisotropies with preferred axes in the universe we exam- 

ine the consequences of the simple hypothesis that the three-dimensional space has a 
^ ! global TZ^ X topology. We take the radius of the compactification to be the ob- 

^ I served cosmological scale beyond which the accelerated expansion starts. We derive 

the induced corrections to the Newton's gravitational potential and we find that for 
distances smaller than the radius the leading 1/r-term is corrected by convergent 
power series of multipole form in the polar angle making explicit the induced anisotropy 
by the compactified third dimension. On the other hand, for distances larger than the 
5^ . compactification scale the asymptotic behavior of the potential exhibits a logarithmic 

dependence with exponentially small corrections. The change of Newton's force from 
to 1/r behavior implies a weakening of the deceleration for the expanding universe. 
Such topologies can also be created locally by standard Newtonian axially symmetric 
mass distributions with periodicity along the symmetry axis. In such cases we can use 
our results to obtain measurable modifications of Newtonian orbits for small distances 
and flat rotation spectra, for large distances at the galactic level. 



1 Introduction 



Nowadays there is accumulating evidence that the standard cosmological model at 
very large distances might not be fully consistent with the astronomical data [T]-|llj. 
Indeed, observations at very large scales point to the existence of a possible anisotropy 
in the global structure of the universe. Similarly, cosmological data for quasars indicate 
that the linear polarization angles of the emitted photons are correlated to a preferred 
axis at huge distances of the order of 1.5 Gpcs [2J. In fact, there is evidence for a 
preferred cosmological axis when comparing accelerations and large peculiar galaxy 
velocities at different directions. Further, from CMB anisotropics observed at the 
lowest multipoles in the spherical harmonic decomposition of temperature, there is 
a hint of the existence of preferred axis[T]. The explanation of rotational curves of 
galaxies implies the possible existence of dark matter but eventually there are other 
explanations [9", T!0]. The standard tests on the validity of Newton's law and general 
relativity gravitational corrections are confined inside the solar system, binary pulsars, 
and gravitational lensing. 

The hypothesis that the Universe has a non-trivial topology and it is possibly multi- 
connected has been considered since the beginning of the 20th century }14|^ . Further, 
similar ideas were developed afterwards |12j where a non-trivial topology was also 
assumed in considerations with regard to the quantum creation of the Universe and the 
origin of inflation [16^ [TT] . 

In the present note, we consider the possibility that matter does not only influence 
the surrounding space-time geometry but also the topology in various scales. This idea 
comes from the well known fact, that in order to describe local physics we need to select 
a complete set of smooth functions on the whole space in order to form delta functions. 
Greens functions etc, but then the topology of the space is crucial in determining this 
complete set. One possible and simple approach is to consider violations of Newton's 
law emerging from local changes of the topology of the Euclidean space TZ^ to TZ^ x S^. 
In a subsequent work we will develop further this idea in the framework of Einstein's 
Theory of General Relativity. At present, we will confine our analysis on the Newtonian 
limit of the above idea i.e. the global topology change of TZ"^ Euclidean space to TZ^ x S^. 
We will obtain then, the exact solution for the Newton's potential interpolating between 
short distances where the usual ^-behavior appears, to the large distance behavior logr 
for distances much larger the radius RofS^. We will see that we obtain axial symmetric 
solutions with multipole terms and indications of correct qualitative behavior for the 
rotational spectra of galaxies. Other tests like gravitational lensing and spectra of 
CMB [7] require further studies. 

We find also that our result is equivalent to that obtained from an infinite discrete 
and periodic mass distribution along the decompactified dimension. An analysis of 
these implications on the Keplerian orbits of Stars at distances where the first non- 
leading corrections are important, is also presented. A detailed comparison with the 
astrophysical data and the observed uncertainties in beyond the scope of this letter 

^For recent surveys see for example [15| . 
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and will be the subject of a future publication. In the present note, in section 2 we 
derive and discuss in detail the properties of the TZ^ x potential and in section 3 we 
analyse its implications on the trajectories of massive objects. The results are discussed 
in section 4. 



2 Gravitational Potential in 'R? x 

Models for modifications of Newton's Gravity have already been proposed for large 
as well as for minuscule distances, the latter being only relevant in the presence of 
extra compact dimensions [21J. A survey of recent possibilities of modifications of 
the Newtonian Gravity can be found in \10\ II ij^ . Here, in this work we perform a 
mathematical analysis which will prove useful to examine whether in the absence of 
Dark Matter, such modifications could also be interpreted as a possible effect of the 
ordinary matter in the topology of the space. 

We assume a three dimensional space with the topology TZ^ x S^. We parametrise 
the compact dimension of 'radius' R as z = Rq with q = [0, 27r). The Poisson equation 
for a point unit mass in this case is written 

Vi+i"^ = 5^^\r-To)^5{q-qo) (1) 

with T = S M^. Expanding in Fourier modes, the solution can be expressed as 

follows 



oo ^ 2 



m=-oci'' 

Performing the k Gaussian integration the potential takes the form 

1 

'Sn'^R 



1 OO OO 2 2 J 



m=— oo 



To proceed with our computation, we split the infinite sum ([3|) into two parts 
1 



Stt^R 



°° ds ^ r°° ds 



T~ ub \ ^ ; _i gin. 

e 4s \-2 cos mq e 4s r 



m=l 



^0 + 6^ (4) 



The first integral diverges, however, we observe that it is the result of the Fourier 
integration of the zero-mode 



^For other types of modified gravity scenarios related to dark matter, the accelerating expansion of 
the universe etc see for example [18. ,19. ,20j . 
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This obeys the two-dimensional Poisson equation 



with the known solution 



2^ 



— i In- 

47r2 R To 



(6) 



(7) 



The result of the second integration in @ is denoted with 5^ and it equals twice 
the modified Bessel function 2KQ{mT/R). Hence, adding the two contributions the 
potential takes the form 



Att'^R 



In 2 cos (m q) {rriT / R) 



(8) 



We can obtain a closed form of the second term using a known formula from [22j . 
Plugging in the appropriate variables, we get 



^ It 

Ko{m-) cos mq= 2 + 4^;^^) + 



ttR 



771=1 



vr \ - 

+ 9E 



1=1 



4vri?' 2Vr2 + (i?g)2 
R 1 



+ 



R 



v/T2+i?2(27r/- g)2 27r/ y/T2 + i?2(2vr/ + g)2 



27r/ 



Using this expansion we notice that the logarithms in 
further choose the constant tq to be tq = 4TTRe~"' ~ 
the potential valid for all distances takes the form 



8|) cancel each other. We may 
2.246 vri? so the final form of 



1 



1 



4vr + R2 q2 



^ 00 

-y 

1=1 



+ 



1 



V'r2+i?2(2vr/- qf InlR + R2(^27rl + qf 



2ttIR 



(9) 



Prom this latter expression, we observe that the three dimensional Newton's potential 
is recovered in the limit R — )• 00, (7—7-0 with z = Rq finite. We further observe that 
form of the potential in ([9]) can also be represented as the Newtonian potential of an 
infinite discrete mass distribution of period 2ttR along the decompactified dimension. 



2.1 The potential at short distances 

The above form of the potential encodes the axial symmetry of the topology of space 
and therefore it is interesting to investigate particular limiting cases with respect to 
the coordinate axes. However it is not easy to elaborate it further in its present form. 
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In order to do this, we need to introduce a convenient parametrization. To this end, 
we define 

r 



and we note that z = Rq = r cos 9. Now we can expand the square roots into Legendre 
polynomials for small r/{2TiR) 

1 1 1 



^r2 + i?2(27rZ± 



2ttIR 



1 + a2 lb 2ai cos ( 



2-kIR 



Y,Pk{±cose)a\ 
k=0 



Substituting into the series form of $ while using the property Pk{— cos ( 
and changing the order of the two summations we obtain 

„2n 



(11) 



(-)^Pfc(cOS^) 



k 1 

4tt r 



k 



$ = __ 



oo 



1 



27r ^ 2-kIR 
1=1 



E 

.n=0 



{2TrlR) 



2n 



P2n(cOSe)-l 



-Ai l + 2j:C(2n + l) 

V n=l 



2n+l 



27rR 



j P2„(cos^) 



(12) 



We remark first that the above expression for the potential holds for small distances 
r/{2'KR) < 1 and moreover care should be taken in the interpretation of the spherical 
coordinate variables r,0 in order to implement the topology of the space Ti? x S^. 

In the following we are going to examine the behavior of the potential in various 
directions of TZ^ x space. We study two cases for the potential, the first along the 
z-axis (9 = 0) and the second on the {x,y) plane (0 = |) . 

In the first case, {9 = 0), we have -Pn(l) = 1 and the infinite sum in (jl2p can 
be expressed in terms of the Euler constant 7 and PolyGamma functions so that the 
potential along the z-axis is given in closed form 



V{r) 



k k 
-+-/ + — 
r ttR 



(-3+*'"'(^n)) 



27rR 



(13) 



For the second case ^ = f and 



P2ni0) 



(-l)"(2n)! 
22«(n!)2 



The potential in this case takes the form 



Vir) 



00 

I + 2E 

n=l 



(-l)"(2n)!C(2n + 1) 



22" (n 



n2 



2n+l 



2ttR 



(14) 



For large n the coefficients of the powers of the expansion parameter s = r ji^ixR) goes 
as (— l)"/T/n, so we can approximate the above sum by appropriate combinations of 

the polylogarithmic functions Lii. The above considerations describe the behavior of 

2 

the potential near the decompactification limit r — )• 0. Figure [T] shows the modification 
of the gravitational potential for 7^} x topology in comparison to TR? case. 



5 



2.2 The potential for large distances 

Because of the compactification radius bound, large distances where r/{2iTR) > 1, can 
be reached only in the direction 9 = ^ or q = 0. In this case, expanding the Bessel 
function in ([8|) for large values of t = mp with p = t/R we have 



Ko{t) 



ra=0 



[(2n-l)!!]2 1 



2t 
vr 



1 



1 32 

+ 



(-8t)" 

(3-5)2 1 



2r V' 8t 2!(8i)2 3! 
and performing the sum in ([8]) we obtain 



, (3-5-7)2 1 



4! 



(8i)^ 



(15) 



[(2n-l)!!]2 (-l)'^ 



m=l 



n=0 



vr 



— Lii (e 
2p V 2 * 



n! 



(Spy 



— Lis (e 
80 2 ^ 



") + — ^Lis (e-") +---^ 



(16) 



The complete potential is obtained by adding the logarithmic term. This in the large 
t/R regime becomes 



Vcip) 



1 1 



2tt^R 

where C is the constant 



C + ilog(p)-yj(Li.(e-)-lL.3(e-) + ...)) 



(17) 



C = ln(2V^) - 7/2 



We note here that the corrections of the logarithmic term are exponentially small. In 
figure [2] we plot the effective potential Ve//(T) = Vc + i'^/{2mT'^), taking three values 
for (./{2m) = 0.06,0.09,0.12 where £ is the conserved angular momentum and m the 
mass. For large values of t/R we add correction terms Li„_|„i up to n = 2. These 
correspond to the right piece of the curves (red color). At distances around r ~ 2ttR 



VxyW 




rlRc 



r/Kc 



Figure 1: Modifications of the Gravitational potential (in arbitrary units) a) along the 
z-axis b) along the xy-plane. Blue (lower) lines corresponds to T^^-case, the red ones 
(upper) to 7^2 x topology. 
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the curves cross the abscissa and for larger r-values, the potential becomes positive and 
increasing. Clearly this is due to the logarithmic dependence of the potential. To plot 
the potential for small r/i? we use the formula (|12|) taking n = 10 terms of the sum 
(in practise only few of them suffice). These values are represented by the left part of 
the curves (blue color). We observe that the two curves match perfectly around the 
area /) ~ 1 , thus our approximation is valid for small r / R < 1 as well as large r/i? > 1 
values. 

We observe now that the corrections terms constitute an alternating series of poly- 
logarithmic functions of the exponential e~P. Thus increasing the values p > 1, their 
value decreases and the dominant contribution comes from the logarithmic term. So, 
for large distances the rotational spectra of massive objects are flat 

V constant 

We notice here that axially symmetric periodic Newtonian mass distributions along the 
^-direction can mimic the effect of the compactification 7i? x 5^. So although presently 
there is no observational basis for such mass distributions, ( for example concentration 
of galaxies along cosmic strings) we can choose the compactification radius R to be of 
the order of the average intergalactic scale. This way, our mechanism could mimic the 
effects of the dark matter. 

3 Modification of the orbits 

The modifications of the potential discussed above, are expected to be negligible at 
small distances. Indeed, one can notice that the relative coefficient of the first correction 
term is ^t3'-P2(0) ~ —0.005 , thus for r ^ i? its contribution diminishes. There are 




Figure 2: Plot of the approximate formulae of 7^^x5^ gravitational potential Vef f{T / R) 
for small t/R { left blue curves) and large t/R (right red curves). 
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considerable deviations from Newton's mechanics however, when deahng with distances 
at very large scales. In this case the relation ([8]) gives the leading logarithmic behavior 
with exponentially small corrections. 

At intermediate distances the effect of the corrections will be important but we need 
to know the exact behavior of the potential which, although complicated algebraically, 
numerically interpolates nicely between the 1/r and log(T) behavior. Here we examine 
the effects of the violations of the Newton's law in the R^-plane for the orbits of massive 
objects and for small distances. 

For distances inside the compactification scale, where the first non-leading cor- 
rection to the Newtonian law plays a role it is interesting to consider the nature of 
modifications on the orbits of massive objects. We start with the calculation of the 
force taking into account the first non- leading term and in polar coordinates (r, 6, (p) 
we find 



F-F - J ^ I 3K(3) 



The form of the correction to Newton's force is similar to the tidal forces on a dipole. 
At this level of approximation we observe that the axially symmetric additional forces 
are small and repulsive for cos(20) > — 5 and attractive otherwise. For the study of 
the effect of this correction to the orbits we restrict our attention to the planar motion 
{6 = f ), Fg = while 

The corresponding equation of motion is 

cfu km f 2C(3) 1 \ 

where u = ^. This equation is easily solved using the conservation of the angular 
momentum, 

V-Vo = / , = (21) 

where 

£ = mr'^ip (22) 

The correction to the Newton's potential is an even-power series of r, the most dominant 
being the harmonic term. Thus, the effective potential can be approximated by 

K„(r) = -'-U-J^-Sm.(LY) (23) 

' r V 2mkr (27r)3 \RJ J ^ ' 
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where 6 = for the classical Newton's law and 5 = 1 for the corrected potential. 
The exact computation of the integral would introduce complicated elliptic functions, 
however, since the correction is small it can be treated perturbatively by expanding 
appropriately the denominator and computing the relevant integrals. Instead of this, 
we will work out in the next section the corresponding equation of motion. 



3.1 The equation of motion 

In the approximation discussed above we may analyse easily the motion using the 
standard approach by deriving the differential equation for the new variable u = 1/r 
introduced above. The perturbed differential equation is written as follows 

^ + u = a{l + eu-^) (24) 

where a = km/i'^ and the parameter e can be deduced from (pOj) 

^_ C(3) 
47r3i?3 

We assume solutions of the form 

u = uo + eui (25) 
where uq is the classical elliptic solution 

, , o(l — e^) 1 , , 

1 + e cos[ip) uo[(pj 

where as usually, a, e stand for the semimajor axis and the eccentricity We seek 
solutions for the perturbed equation which have the same angular momentum as the 
unperturbed one, thus Si = 0. Substituting into (p^ and using the equation (p6]) for 
Wo ((/?), one finds 

This has the solution 

Mv) = = 2 2^ 

where Isi'^) and Id'^) are given by 

, 3e 1 / 11- e ip\ (2 + 6^ + (l + 2e2) ecosv?) sin(/? 

Urn) = FTTTtan"^ \ tan- - ^ — . ^ ^ 

{1-e^f^ yll + e 2 J 2(l-e2)2(l + ecos<^)2 

Isif) = — 



2e V(l + e)2 (1 + ecosc^)^ 



''These are defined as a = — and e — y I ~ ~ v ^ ~ Tcma' ^'^^^ ^1^° that the parameter a 
is related to the physical parameters a,e by = a(l — e^). 
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Figure 3: The ratio ^^^^ with r being the distorted orbit in the presence of the 
correction term in TZ^xS^ topology. The three curves (upper, middle, lower) correspond 
to three different ratios a/R = 1/6,1/5,1/4 where a is the semimajor axis of the 
unperturbed orbit and R the 5^-radius. 

To get a feeling of the corrections implied on the trajectories of physical objects we form 
the ratio of the corrected over the unperturbed solution. To this end, we expand the 
ratio in terms of the dimensionless quantity e/a^ and find that for small eccentricities 
e <^ 1 the essential modification of the orbits is captured by the formula 



For R of the order of galactic distances and a a typical solar system distance, this 
implies a negligible distortion on the orbits. 

Next we study the impact on the trajectories of massive objects with bounded 
motion, using three values for the ratio A = 1/6, 1/5, 1/4 while assuming the eccentricity 
e = 0.068. For an assumed massive object subject to the forces due to a potential 
generated at the origin, the corresponding ratios are plotted in figure ([3]). It is obvious 
that for massive objects at the boundaries of galaxies the radius R has to be taken of 
the same order to have an observable modification of Newtonian orbits. 

We close our analysis with the following remark. Suppose that in the absence of 
the TZ'^ X 5^-effects the orbit of an object with eccentricity e is determined as usual by 
ro{ip) given by (f26l) . Substituting, the solution ([291) can be also written as follows 




(29) 



where ro{ip) is given by ([26j) and we have defined 



/c(3)y _a_ 



\ 2 J ttR 




1 - e sin2 f 
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For objects satisfying = 2e/(l + e), or equivalently 



a = TT R 



( 



C(3)l + e 



1 4e 



) 



3 



and to first order approximation on the eccentricity e, we obtain 



r{(/)) = a(l — e) = constant 



Therefore, certain "would be" elliptic orbits in TZ^ space, in the case of TZ^xS^ topology 
appear to be approximately circular. 

4 Discussion and conclusions 

Motivated by the fact that several recent astronomical and cosmological observations 
are challenging the theory of General Relativity at very large distances, due to accel- 
erated expansion, axial anisotropics, planarity of CMB multipoles, dark energy etc in 
this note we have considered the role of a non-trivial topology of the universe in the 
Newtonian approximation without invoking extra dimensions. 

Assuming the experimentally determined flatness of the spatial universe, we test 
at a qualitative level the proposal of the simplest periodic multiverse along the third 
spatial dimension or equivalently the assumption of a spatial topology TZ^ x with a 
radius of the compactified dimension of the order of the distance where the acceleration 
starts, L = 2itR. 

Within the aforementioned topology of the space, we have determined the exact 
Newtonian potential at any distance r. 

For small r/R ratios we find the corrections to the three dimensional Newton's 1/r- 
potential to be of the form of tidal forces (r/i?)" P„(cos^), with an anisotropy axis 
along the third dimension. Another interesting finding at these distances is the repulsive 
accelerated motion along the axial direction in consistency with the observed large 
peculiar velocities of galaxies along the observed axes at some statistical significance 
level. 

On the other hand, for large distances, where r/R> 1, we get a smooth transition of 
the potential to the two dimensional one, log(r), with exponentially small corrections. 
This has the immediate consequence that for distances larger than L we obtain smaller 
decelerations due to the appearance of the two dimensional attractive force 1/r, instead 
of the 1/r^ of three dimensions, thus a relative acceleration. 

Moreover, we expect that for such large distances or angular separations there must 
be a flattening of the CMB multipoles due to lower dimensionality of space. 

Choosing smaller scales for the compactification radius could be afforded only in 
some General Relativity model, where the compactification is local and not global. 
Another possibility is to mimic the compactification of the third dimension with spatial 
axially symmetric and periodic in the third direction Newtonian mass distributions. In 
such a patch -say around a galaxy- we obtain that the trajectories of massive stars 
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at the boundaries of the galaxy have flat rotational spectra and small but measurable 
perturbations of the periodic orbits for smaller distances. 

We find the above results encouraging and the proposed topological modification of 
gravity worth of extending, in a future work, in the framework of General Relativity. 
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